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1.  INTRODUCTION  AND  BACKGROUND 


The  objective  of  this  program  is  to  develop  both  a  qualitative  and 
a  quantitative  understanding  of  the  fluid  mechanics  and  heat  transfer 
mechanisms  that  underlie  the  measured  performance  of  the  spirally 
fluted  tubes  under  development  at  General  Atomic.  The  reason  for  the 
interest  in  the  spirally  fluted  tubes  is  that  results  to  date  have 
indicated  three  advantages  to  this  tubing  concept: 

1.  The  fabrication  technique  of  rolling  flutes  on  strip  and 
subsequently  spiralling  and  simultaneously  welding  the  strip 
to  form  tubing  results  in  low  fabrication  costs* 
approximately  equal  to  those  of  commercially  welded  tubing. 

2.  The  heat  transfer  coefficient  is  increased  without  a  concomi¬ 
tant  increase  of  the  friction  coefficient  on  the  inside  of 
the  tube.  In  single-phase  axial  flow  of  water,  the  helical 
flutes  continuously  induce  rotation  of  the  flow  both  within 
and  without  the  tube  as  a  result  of  the  effect  of  curvature. 
This  rotation  enhances  turbulent  exchange  both  on  the  inside 
and  outside  of  the  tube  in  the  immediate  vicinity  of  the 
wall.  When  heat  is  transferred  from  a  high- temperature  fluid 
external  to  the  tube  through  the  tube  wall  to  a  cold  fluid 
within  the  tube,  the  density  gradients  are  such  as  to  enhance 
turbulent  exchange  in  both  flows.  The  improvement  in  the 
heat  transfer  coefficient  results  from  the  enhancement  of  the 
turbulent  exchange. 

An  additional  effect  from  rotation  of  the  fluid  inside  of  the 
tube  has  been  observed  for  a  helix  angle  of  30°.  The 
friction  coefficient  does  not  increase  with  the  increased 
heat  transfer;  in  fact,  it  is  somewhat  lower  than  that  of  a 
smooth  round  tube  at  equivalent  Reynolds  number. 
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The  combination  of  these  two  effects  results  in  the  achieve¬ 
ment  of  values  considerably  in  excess  of  one  for  the  ratio  of 
the  turbulent  exchange  coefficients  of  heat  and  momentum  (the 
reciprocal  of  the  turbulent  Prandtl  number).  This  ratio  is 
equivalent  to  the  ratio  of  twice  the  Stanton  number*  multi¬ 
plied  by  the  Prandtl  number  raised  to  the  two-thirds  power, 
to  the  friction  coefficient. 


3.  An  increase  in  condensation  heat  transfer  on  the  outside  of 
the  tube  is  achieved.  In  a  vertical  orientation  with  fluid 
condensing  on  the  outside  of  the  helically  fluted  tube,  the 
flutes  provide  a  channel  for  draining  the  condensed  fluid. 

The  surface  tension  forces  draw  the  condensate  film  from  the 
crests  of  the  flutes  into  the  troughs.  This  results  in  the 
major  portion  of  the  crest  having  a  very  thin  film*  which 
greatly  reduces  the  resistance  to  heat  flow  through  the  crest 
area.  The  result  of  providing  the  drainage  in  the  troughs 
and  decreasing  the  resistance  to  heat  flow  over  the  crests  is 
a  substantial  improvement  in  heat  transfer  performance  over 
conventional  smooth  tubes.  This  effect,  which  was  first 
proposed  by  Gregorig  [1].  is  not  novel  to  this  tubing 
concept.  In  fact  it  is  the  basis  for  a  number  of  different 
proposed  enhanced  tubes,  some  of  which  have  straight  flutes 
and  some  spiralled.  However,  the  combination  of  this 
enhancement  on  the  outside  of  the  tube  in  condensation  with 
the  enhancement  on  the  inside  of  the  tube  in  single-phase 
water  flow  can  lead  to  significant  improvements  in  heat 
transfer  performance. 

The  manufacturing  development  of  the  spirally  fluted  tube  Is  underway  and 
its  evaluation  will  be  reported  on  later.  Th?  basis  of  the  other  two 
advantages  follow: 
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Qi.  Helically  Fluted  Tubes  in 


Water  Flow 


The  heat  transfer  performance  on  the  inside  of  helically  fluted 
tubes  as  indicated  by  the  Nusselt  modulus,  Nu/PrO*^,  as  a  function  of 
the  Reynolds  number  is  shown  in  Fig.  1.  The  Nusselt  number  is  defined 
as  Nu  =  hdjjyj/k,  where  h  is  the  conductance  coefficient 
(Btu/hr-ft2-°F) ,  d^yd  is  the  hydraulic  diameter  (ft.)  and  is  equal  to  4 
x  cross-sectional  area  (f t2)/perimeter  (ft.),  and  k  is  the  thermal 
conductivity  of  the  water  (Btu/hr-°F-ft. ) .  The  Reynolds  number  is 
defined  as  vdhy^/v,  where  v  is  the  mean  axial  velocity  (ft/sec)  and  v 
is  the  kinematic  viscosity  (ft2/sec).  The  Prandtl  number  is  Pr  =  c  p/k, 
where  cp  is  the  specific  heat  at  constant  pressure  and  p  is  the 
absolute  viscosity.  The  Dittus  Boelter  correlation,  Nu/Pr®*^  =  0.023 
Re®*®,  for  a  smooth  tube  is  shown  for  comparison. 

The  ratio  of  the  ordinates  for  the  four  curves  shown  in  Fig.  1  to 
that  of  the  Bmooth  tube  curve  indicates  the  amount  of  heat  transfer 
enhancement  achieved  on  the  inside  of  the  helically  fluted  tubes 
relative  to  a  smooth  tube.  These  curves  do  not  include  the  effect  of 
the  extended  area  of  the  helically  fluted  tubes  relative  to  that  of  a 
smooth  tube  of  the  same  cross- sectional  area.  This  is  an  additional 
factor  of  1.33  for  the  tube  reported  here;  in  some  of  the  tubes  under 
development  the  extended  area  is  1.6  greater  than  the  area  of  a  smooth 
tube  of  equal  cross-sectional  area.  Curves  1,  3,  and  4  are  the  results 
from  tests  conducted  by  Heat  Transfer  Research,  Inc.  (HTRI),  on  a 
l-in.-o.d.  tube  that  had  flutes  with  a  30°  helix  angle  and  was 
approximately  7  ft.  long.  (The  helix  angle  is  defined  as  the  angle  the 
helix  makes  with  the  tube  axis.)  Curve  1  is  for  cooling;  that  is,  the 
direction  of  heat  flow  is  radially  outward.  Hot  water  flowing  on  the 
inside  of  the  tube  was  cooled  by  colder  water  flowing  on  the  outside. 
Although  enhancement  is  achieved  in  this  mode  of  operation,  it  is  less 
than  in  the  other  curves  since  the  density  gradient  is  stabilizing. 


3 


NUSSELT 

MODULUS 


H 

I 

I 


Curves  3  and  4,  however*  show  significant  enhancement.  These  data 
were  obtained  with  the  direction  of  the  heat  flow  radially  inward* 
which  leads  to  density  gradients  that  are  destabilizing  and  hence 
increases  turbulent  exchange.  Steam  was  condensed  on  the  outside  of 
the  tube  by  colder  water  flowing  on  the  inside.  The  condensing 
pressure  was  approximately  9  atm  for  curve  4  and  5  atm  for  curve  3  in 
these  tests.  The  steam  flow  was  parallel  to  the  tube  axis*  which  was 
horizontal.  Data  from  two  tests  conducted  by  D.  J.  Reilly  at  the  Naval 
Postgraduate  School  [2]  on  a  5/8-in.  30°  helix  angle  spirally  fluted 
tube  are  shown  by  curve  2.  In  these  tests  steam  at  0.25  atm  flowed 
transverse  to  the  tube  axis,  which  was  horizontal.  Curves  2*  3*  and  4 
indicate  that  as  the  condensing  pressure  and  hence  temperature  were 
increased  relative  to  the  tube-side  coolant,  heat  transfer  increased. 
The  higher  condensing  temperatures  indicate  higher  heat  fluxes  and 
result  in  larger  temperature  gradients  on  the  inside  of  the  tube. 

The  frictional  performance  of  the  30°  helix  angle  tube  can  be  seen 
in  Fig.  2.  where  the  isothermal  friction  coefficient  (Moody  definition) 
is  plotted  against  Reynolds  number  for  both  the  HTRI  tests  (indicated 
by  the  line  and  the  Reilly  tests  (indicated  by  the  triangles).  The 
hydraulic  diameter  is  used  in  the  calculation  for  the  frictional 
coefficient. 

Since  a  decrease  in  the  frictional  coefficient  in  a  tube  in  which 
enhancement  of  the  heat  transfer  is  achieved  is  surprising,  a  series 
of  additional  tests  were  undertaken.  These  experiments  were  designed 
to  eliminate  the  uncertainties  in  the  measurement  of  frictional 
coefficients  in  test  rigs  that  were  designed  for  heat  transfer 
measurements.  The  preliminary  data  from  the  tests  underway  at  the 
present  time  confirm  the  frictional  coefficients  depicted  in  Fig.  2 
over  a  Reynolds  number  range  between  40*000  and  80*000  with  air  as  the 
fluid.  Tests  with  water  as  the  working  fluid  are  scheduled  which  will 
also  extend  the  Reynolds  number  range  of  the  data.  The  croas-sectional 
area  of  the  tubes  varied  from  0.127  in.^  in  the  Reilly  tests  to  0.885 
in. ^  in  the  tests  underway  at  present.  The  cross-sectional  area  of  the 
tube  in  the  HTRI  tests  was  0.576  in.^. 
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1.2  Combined  Heat  Transfer  and  Friction 

While  the  analogy  between  momentum  and  heat  transfer  suggested  by 
Reynolds  and  modified  by  Colburn  is  valid  only  under  rather  restrictive 
conditions,  it  is  almost  by  definition  an  index  of  the  degree  of  heat 
transfer  enhancement  achieved  relative  to  the  frictional  penalty.  It 
is  directly  related  to  the  ratio  of  the  exchange  coefficients  of  heat 
to  momentum.  Therefore,  it  is  of  interest  to  examine  the  heat  transfer 
and  friction  data  in  the  context  of  the  Colburn  analogy.  The  Colburn 
analogy  is  defined  as  2j/f*  and  is  equal  to  unity  for  a  straight  smooth 

tube,  where  j  =  N8Pr2/3.  it  shown  in  Fig.  3  for  the  high  flux  case, 

the  low  flux  case,  the  cooling  case  and  for  Reilly  tests. 

There  are  two  items  of  interest.  First,  and  least  significant,  is 

the  positive  slope  that  results  from  the  negative  slope  of  the  friction 

coeff icient/Reynolds  number  characteristics  of  Fig.  1,  which  is  greater 
than  that  of  a  smooth  tube.  More  important  is  the  level  of  the  Colburn 
factor,  2j/f,  which  achieves  a  value  of  2.3  at  low  flux  and  a  value  of 
2.75  at  high  flux,when  a  value  of  one  is  the  case  for  a  smooth  tube  and 
less  than  one  for  most  other  enhancement  techniques. 

1 .3  Coadea&ation  ol  Helically  Elated  Tubes  la  A  Yextical  Configuration 

Two  tubes  with  different  helical  angles  (one  a  30°  helix  and  the 
other  a  45°  helix)  were  tested  at  Oak  Ridge  National  Laboratory  with 
Refrigerant  11  condensing  on  the  outside  of  the  tube,  which  was  cooled 
by  water  on  the  inside.  The  test  data  depicting  the  performance  of 
these  tubes  are  shown  in  Figs.  4  and  5.  In  Fig.  4,  the  condensing  film 
coefficient  is  the  ordinate  and  the  heat  flux  is  the  abscissa.  The 
increase  in  performance  over  the  smooth  tube  decreases  with  heat  flux 


The  frictional  coefficients,  f*  and  f,  used  in  the  text  and  figures 
are  respectively  based  on  the  Moody  and  Fanning  definitions,  which 
differ  by  a  factor  of  4.  The  Moody  coefficient  is  larger  than  the 
Fanning  coefficient  by  this  factor. 
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45°  HELIX 


=  Q/A,  HEAT  FLUX  (BTU/HR-FT*) 


but  remains  quite  substantial.  Figure  5  shows  the  same  data  with  the 
heat  load  as  the  ordinate  and  the  condensing  temperature  difference  as 
the  abscissa.  Again,  the  gain  over  the  smooth  is  larger  at  the  lower 
heat  loads.  The  45°  helix  appears  to  have  better  performance  at  the 
lower  heat  flux.  It  should  be  pointed  out  that  the  45°  tube  had  fewer 
flutes  and  hence  troughs,  so  the  troughs  would  fill  sooner.  The  liquid 
capacity  of  each  of  the  flutes  in  both  of  these  tubes  was  approximately 
the  same.  Therefore,  a  larger  number  of  flutes  would  drain  more 
liquid.  Another  fact  that  should  be  kept  in  mind  is  that  the  latent 
heat  of  Refrigerant  11  is  quite  low,  so  a  considerably  larger  amount  of 
liquid  must  be  drained  from  the  tubes  than  would  be  the  case  for  steam 
or  ammonia. 

These  data  indicate  that  this  tubing  concept  is  uniquely  suited  to 
condensers,  as  enhancement  is  achieved  in  a  vertical  configuration  on 
both  the  outside  condensing  surface  and  on  the  inside  water  surface 
without  an  increase  in  friction.  These  advantages  can  be  particularly 
significant  when  the  weight  and  volume  of  the  condenser  is  important  as 
it  is  in  naval  vessels  either  for  the  main  propulsion  system  or  for 
bottoming  cycles. 

Any  rational  explanation  of  the  heat  transfer  enhancement  that 
occurs  without  a  frictional  penalty  must  focus  on  the  fundamental  fluid 
mechanical  processes  of  turbulence.  It  is  widely  recognized  that  the 
resistance  to  the  transmission  of  heat  (or  any  other  scaler)  from  the 
wall  to  the  internal  flow  in  a  tube  is  in  the  laminar  sublayer  adjacent 
to  the  wall  where  the  process  is  dependent  on  the  molecular  conductivity. 
Improvements  in  the  rate  of  radial  heat  transport  resulting  from 
increases  in  the  turbulence  level  in  the  flow  normally  results  in 
greater  Increases  in  the  momentum  losses  relative  to  the  increase  in 
heat  transfer.  However,  there  are  examples  of  instabilities  in  the 
atmospheric  boundary  layer  that  result  in  much  greater  heat  exchange 
than  momentum  exchange.  A  laboratory  experiment  that  illustrates  the 
atmospheric  phenomena  is  that  of  Mizushina,  et  al.[3]  The  exchange 
coefficients  of  heat  and  momentum  in  thermally  stratified  flow  in  an 


CONDENSING  TEMPERATURE  DIFFERENCE  AT  (°F) 


open  channel  under  both  stable  and  unstable  conditions  were  determined 
from  measurements  of  the  temperature  and  velocity  gradients.  Their 
results  indicate  that  the  dif fusivities  of  heat  and  momentum  decrease 
with  stability  and  the  ratio  of  the  diffusivity  of  heat  to  that  of 
momentum  decreases  with  increasing  Richardson  number.  Under  unstable 
conditions.  the  dif fusivities  of  both  beat  and  momentum  increase  with 
instability  and  the  ratio  of  the  dif  fusivities  of  heat  to  that  of 
momentum  increses  to  a  maximum  value  of  3  as  the  stratification  shifts 
from  neutral  to  a  weakly  unstable  condition.  It  then  decreases 
gradually  as  the  instability  increases.  It  is  our  view  that  these 
effects  of  buoyancy  in  a  gravitational  field  are  replicated  in  the 
spirally  fluted  tube  by  the  rotational  flow  induced  by  the  flutes. 
There  is  an  axisymmetrical  core  region,  stabilized  by  swirl,  surrounded 
by  an  annular  region  which  is  close  to  the  flutes,  highly  three 
dimensional  and  unstable. 

The  purpose  of  this  research  program  is  to  develop  a  theory  to 
explain  the  increase  in  heat  transfer  relative  to  friction  achieved  in 
the  reported  measurements.  A  combined  experimental  and  theoretical 
program  has  been  underway  at  the  University  of  California  at  San  Diego 
and  the  General  Atomic  Company  and  is  being  expanded  to  include  work  at 
the  University  of  Manchester  and  Imperial  College.  The  initial  work 
is  reported  herein  and  includes  the  construction  of  the  experimental 
facilities,  the  development  of  a  casting  technique  to  provide  a 
transparent  test  section  which  has  an  internal  sprially  fluted  flow 
passage,  the  flow  visualization  studies  performed  to  date,  an 
experiment  to  measure  the  overall  heat  transfer  coefficient  using  air 
as  the  cooling  medium,  the  results  to-date  of  the  theoretical  analysis 
and  the  work  planned  for  the  future. 
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2.  EXPERIMENTAL  PROGRAM 


The  foci  of  the  experimental  research  effort  during  the  past  year 
have  been  the  design  and  construction  of  both  a  water  and  an  air  flow 
facility.  The  water  flow  facility  has  been  used  primarily  to  study  the 
flow  field  in  a  spirally  fluted  tube  by  means  of  flow  visualization 
techniques,  while  the  air  flow  facility  has  been  used  to  determine  the 
average  or  bulk  heat  transfer  rates  for  air  flowing  inside  the  tube 
with  the  heat  flowing  radially  inward.  In  the  next  sections  we 

describe  the  two  facilities  as  well  as  the  major  technological 
accomplishments  in  the  development  and  construction  of  the  two 
facilities.  We  next  discuss  the  experimental  results  and  the  plans  for 
future  work. 


2.1  Water  Flow  Facility 

The  water  flow  facility  has  been  designed  to  provide  a  flow  where 
the  flow  field  in  a  spirally  fluted  tube  can  be  studied  both 
qualitatively  and  quantitatively.  The  facility  is  shown  schematically 
in  Fig.  6  and  consists  basically  of  a  pump,  damping  tank,  flow 
straightener,  test  section,  and  return  line.  The  maximum  flow  rate 
corresponds  to  a  Reynolds  number  in  excess  of  10^  for  a  tube  with  a 
hydraulic  diameter  of  0.6-inch. 

Since  the  flow  facility  is  to  be  used  initially  for  flow 
visualization  studies,  the  test  section  has  to  be  optically 
transparent,  the  inner  surface  must  be  goemetrically  similar  to  the 
spiral  fluted  tubes  which  exhibit  enhanced  heat  transfer,  and  the  outer 
surface  must  be  flat.  The  construction  of  such  a  test  section  is  not  a 
straightforward  task  and  therefore  has  consumed  a  rather  large  percentage 
of  our  effort. 
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2.2  Test  Section 


The  transparent  test  section  is  constructed  by  placing  an  aluminum 
spiral  fluted  tube  in  the  center  of  a  2.25  inch  square  mold.  The  sides 
of  the  mold  are  formed  of  1/A  in.  thick  plate  glass.  The  inside 
surface  of  the  glass  is  covered  with  mylar  film  which  is  attached  to 
the  glass  by  means  of  double  sided  masking  tape.  A  small.  1/8  inch 
slot,  is  left  at  the  corner  two  of  the  glass  plates  and  the  other 
corners  are  sealed  with  RTV  rubber.  It  is  through  this  slot  that  the 
resin  is  poured  into  the  mold.  The  entire  mold  is  filled  with 
Castolite  AP,  acrylic-polyester  resin.  The  hardening  of  the  resin  is 
slightly  exothermic  and  is  accompanied  by  some  shrinkage  -  about  1/8" 
for  each  inch  of  thickness.  Therefore  it  is  important  that  the  curing 
or  hardening  of  the  resin  take  place  slowly  and  uniformly  so  as  to  lead 
to  a  minimum  level  of  internal  stress.  This  curing  schedule  can  be 
accomplished  by  thorough  mixing  of  the  recommended  minimum  amount  of 
hardener  in  the  casting  resin.  Since  the  resin  contracts  during  the 
hardening  process,  it  must  also  be  free  to  move  relative  to  the 
boundary  of  the  mold  in  order  to  reduce  the  magnitude  of  surface 
stresses  and  consequent  surface  cracking.  It  is  for  this  reason  that 
mylar  which  does  not  adhere  to  the  resin  is  attached  to  the  interior 
surface  of  the  glass  boundaries  of  mold.  An  alternate  approach  was 
developed  in  which  the  glass  surface  on  the  interior  of  the  mold  is 
coated  with  a  light  layer  of  hard  wax  but  this  technique  yields  a 
surface  which  is  less  transparent. 

The  shrinkage  of  the  resin  does  lead  to  surface  undulations  with  a 
peak  to  peak  amplitude  of  1/32  to  1/16  of  an  inch.  These  undulations 
lead  to  optical  distortion  and  must  be  removed.  The  surface 
undulations  are  too  large  to  be  easily  removed  mechanically.  For  this 
reason,  the  hardened  casting  is  placed  in  a  second  glass  mold.  The 
interior  of  the  glass  walls  are  covered  with  mylar  as  in  the  previous 
step  and  are  placed  about  1/8  inch  from  the  surface  of  the  casting. 
This  thin  layer  is  filled  with  casting  resin.  Shrinkage  of  this  thin 
layer  also  occurs  during  hardening  and  leads  to  surface  undulations 

with  a  peak  to  peak  amplitude  of  1/6A  inch.  These  undulations  can  be 
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removed  by  wet  sanding  and  then  buffing  with  a  polishing  compound 


The  exterior  surface  of  the  spiral  fluted  tube  which  is  placed  in 
the  center  of  the  mold  is  also  treated  so  that  the  fluted,  interior 
surface  of  the  mold  will  be  smooth  and  free  of  discoloration.  This 
treatment  consists  first  of  machining  the  weld  so  that  the  weld  becomes 
geometrically  similar  to  a  flute,  second,  of  cleaning  the  exterior 
surface  of  the  tube  with  a  weak  solution  of  sodium  hydroxide,  and, 
third,  of  buffing  the  tube  surface. 

After  the  resin  from  both  pours  has  hardened,  the  glass  is  removed 
and  the  casting  is  immersed  in  a  30Z  solution  of  sodium  hydroxide  which 
is  used  to  dissolve  the  aluminum  tube.  This  process  takes  2-3  days. 
The  removal  of  the  aluminum  could  be  accomplished  in  a  shorter  period 
of  fime  through  the  use  of  a  more  concentrated  and/or  heated  solution. 
However,  the  heat  release  would  be  higher  and  would  lead  to  a  higher 
temperature  of  the  casting  resulting  in  discoloration. 

Teat  sections  up  to  eight  feet  in  length  have  been  cast.  These 
long  castings  are  mechanically  strengthened  by  placing  1/4  inch 
diameter  aluminum  rods  at  the  four  corners  of  the  mold  which  are 
fastened  to  plexiglas  flanges  at  the  two  ends  of  the  molds. 

The  test  sections  produced  using  the  process  described  in  the 
preceding  paragraphs  are  of  high  optical  quality  which  is  required  for 
the  flow  visualization  studies  and  the  possible  quantitative 
measurement  of  the  flow  field  by  means  of  a  laser  doppler  velocimeter. 

2.3  Flow  Visualization  Studies 

Flow  visualization  studies  were  made  with  the  water  tunnel  shown 
in  Fig.  6  with  a  four  foot  long  test  section  that  was  produced  using 
the  casting  technique  previously  described.  Two  plane  mirrors  were 
attached  to  the  test  section  at  approximately  a  43°  angle.  Figure  7 
shows  the  casting  with  the  two  plane  mirrors.  This  enabled  the 
observation  of  a  top,  bottom  and  side  view  from  one  side  when  the 


lighting  source  was  from  the  opposite  side.  This  turned  out  to  be 
unsatisfactory  as  the  light  intensity  was  much  greater  in  the  direct 
view  (side)  than  in  the  reflected  views  (top  and  bottom).  The  results 
were  improved  by  introducing  light  through  the  back  side  of  the  test 
section  and  through  one  of  the  mirrors,  the  bottom  one.  and  observing 
the  top  and  side  of  the  test  section.  Observations  were  also  made 
through  a  transparent  plate  at  the  end  of  the  water  tunnel  with 
illumination  from  the  plenum  section,  upstream  of  the  test  section. 
Colored  dyes,  red  and  green,  were  introduced  into  the  test  section  at 
five  radial  stations;  wall  (within  the  flutes),  at  the  top  of  the 
flutes,  at  1/8  diameter  at  1/4  dimeter,  and  in  the  center  of  the 
spirally  fluted  tube,  1/2  diameter.  Both  high  speed  motion  and  still 
photographs  were  taken.  Figure  8  is  a  still  photograph  taken  from  the 
side  of  the  test  section.  Red  dyt  is  injected  at  the  top  of  the  test 
section  and  green  dye  at  the  bottom  with  the  water  flow  from  the  right 
to  the  left.  The  dye  injection  was  continuous.  The  flute  spiral  is 
left  handed  and  is  curved  f.om  the  bottom  to  the  top  and  the  far  side 
of  the  test  section  can  be  seen  as  light  hatching  normal  to  the  spiral 
on  the  viewers  side.  The  rotational  nature  of  the  flow  can  be  seen  as 
the  green  dye  rotates  from  its  point  of  injection  at  the  bottom  toward 
the  viewer  and  curves  around  the  tube  while  the  red  dye,  injected  at 
the  top.  curves  away  from  the  viewer  to  reappear  at  the  bottom 
approximately  one  pitch  distance  downstream.  The  two  fluids  appear  to 
be  completely  mixed  in  two  pitch  lengths.  The  spiral  angle  is  30°.  the 
tube  diameter  approximately  1  inch  and  the  pitch  approximately  3 
inches. 

End  views  of  the  flow  in  a  spirally  fluted  tube  can  be  seen  in 
Figures  9  and  10.  These  are  collages  of  individual  photographs  taken 
through  the  transparent  downstream  tunnel  wall  looking  upstream  into 
the  spirally  fluted  tube.  Lighting  was  provided  by  reflecting  light 
from  a  mirror  mounted  in  the  upstream  plenum  section.  Colored  dye  was 
injected  into  the  tube  through  hypodermic  tubes  mounted  on  the  sides 
of  the  transparent  test  section.  A  timing  system  was  designed  to 
inject  the  dye  and  then  after  a  predetermined  delay  release  the  camera 

shutter.  The  system  can  accommodate  two  cameras,  one  viewing  from  the 


Flow  visualization  -  side  viev.  Red  dye  injected  at  top, 
n  dve  at  bettor.,  water  flow  is  from  left  to  right. 


Fig,  10  F'Jnv  visualization,  end  view. 


end  and  the  other  from  the  side.  A  steady  light  source  was  used. 
Individual  photographs  were  taken  to  simulate  the  fluid  motion 
sequentially.  Each  photograph  was  taken  with  an  individual  injection 
of  dye  with  a  different  time  delay.  The  results  taken  to  date  are 
qualitative.  In  Figure  9  the  injection  of  red  dye  was  at  the  flute  - 
core  boundary  and  the  green  dye  at  the  center  line  of  the  tube.  Note 
the  peripheral  spread  of  the  red  dye  as  compared  with  the  green  dye 
particularly  in  picture  1  where  the  injections  can  be  seen.  (Rotation 
is  clockwise.)  In  Figure  10  red  dye  injection  is  at  1/4  of  the  tube 
diameter  and  the  green  dye  on  the  center  line.  In  picture  1  of  Figure 
10  note  the  diffusion  and  rotation  from  the  injectors.  The  rotation  is 
not  as  marked  as  in  the  previous  figure  with  red  dye  injection  at  the 
flute  -  core  flow  boundary.  While  these  photographs,  other  still 
photographs*  the  high  speed  motion  pictures  and  visual  observations 
qualitatively  support  the  notion  of  the  flow  rotation  and  turbulent 
exchange  between  the  flute  and  core  flows*  they  did  not  lead 
to  definitive  quantitative  information.  Since  a  steady  light  source  of 
limited  intensity  was  used,  the  exposure  time  and  resulting  resolution 
were  limited.  This  proved  to  be  a  more  severe  handicap  in  the  motion 
pictures.  Furthermore,  the  inertia  in  the  dye  injection  system  leads 
to  errors  when  sequential  still  photographs  are  compared.  We  anticipate 
that  these  problems  will  be  resolved  with  the  completion  of  the  high 
intensity  xenon  flash  system  that  is  under  development.  This  will  enable 
high  speed  motion  pictures  with  frame  by  frame  flash  synchronization  with 
exposure  times  of  less  than  1/2  of  a  millisecond. 


2.4  Quafttitatiye  Flow  Measurements 

A  necessary  synergism  exists  between  the  quantitative  measurement 
of  velocity  and  the  theoretical  model  which  will  be  used  to  predict  the 
flow  field.  Thus*  measurements  of  the  velocity  field  in  the  spiral 
fluted  tube  are  not  only  of  fundamental  interest  but  are  also  required 
to  provide  insight  and  guidance  to  the  development  of  a  valid 
theoretical  model  of  the  flow. 
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Two  means  of  measuring  the  velocity  are  of  interest:  constant 
temperature  anemometry  and  laser  doppler  anemometry.  The  latter  offers 
several  advantages  in  comparison  to  the  former;  the  spatial  resolution 
is  better  and  the  measurements  obtained  are  neither  contaminated  by 
temperature  fluctuations,  nor  by  gas  bubbles.  Unfortunately,  the  cost 
of  a  laser  system  is  about  ten  times  that  of  the  corresponding  constant 
temperature  anemometer  system.  The  constant  temperature  anemometer 
system  with  spatial  resolution  of  1  in  is  adequate  in  all  regions  of 
the  tube  except  near  the  flutes.  There,  a  laser  system  which  can 
resolve  scales  to  0.2  mm,  would  have  to  be  used.  If  measurements  in 
the  flute  region  become  critical  to  the  program,  purchase  of  a  laser 
system  will  be  proposed.  However,  during  the  next  phase  of  the 
research,  we  propose  to  make  quantitative  flow  measurements  by  means  of 
a  hot  film,  constant  temperature  anemometry  system.  We  next  consider 
some  of  the  problems  associated  with  the  use  of  the  constant 
temperature  anemometry  techniques  in  water  and  then  discuss  the 
proposed  measurements. 

2.4.1  Measurement  Problems  in  Water 

The  technique  is  based  on  the  physical  principal  that  the  rate  of 
heat  transfer  from  a  small  heated  cylinder  is  uniquely  related  to  the 
velocity  of  a  constant  property  fluid.  In  constant  temperature 
anemometry  the  temperature  of  the  heated  cylinder,  which  consists  of 
a  thin  layer  of  platinum  which  is  sputtered  on  a  25y  diameter  by  1  mm 
long  quartz  fiber,  is  held  constant  by  appropriate  electronic  circuitry. 
The  power  required  to  hold  the  temperature  constant  is  measured.  The 
power  depends  on  the  velocity  and  on  the  surface  conditions  at  the 
water-sensor  interface  which  are  unique  to  each  sensor.  Therefore  each 
sensor  must  be  individually  calibrated  -  nominally  on  a  daily  basis  - 
in  a  calibration  tunnel  where  various  stable  flow  velocities  can  be 
obtained.  Contrast  this  with  the  laser  doppler  anemometer  which 
maintains  its  calibration  independent  of  the  properties  of  the  fluid. 

A  second  problem  encountered  in  hot  film,  constant  temperature 
anemometry  is  outgassing  of  dissolved  gas  onto  the  surface  of  the  film. 
This  outgassing  is  due  primarily  to  the  facts  that  the  sensor  is 


heated,  that  the  surface  serves  as  a  nucleation  site  for  dissolved  gas, 
and  that  a  relatively  low  pressure  is  created  as  the  fluid  moves  around 
the  cylindrical  surface  of  the  sensor.  The  occurrence  of  gas  bubbles 
on  the  surface  of  the  film  dramatically  effects  the  heat  transfer  rate 
and  power  required  by  the  sensors.  Thus  great  care  must  be  exercised 
in  removing  dissolved  gases  from  the  water  and  removing  those  gas 
bubbles  that  collect  at  the  sensor  surface.  Much  of  the  dissolved  gas 
can  be  removed  by  passing  the  water  through  a  nozzle  followed  by  a 
sudden  expansion  where  suction  is  created.  Gas  bubbles  can  be  blown 
off  the  sensor  surface  by  means  of  a  water  or  air  jet  and  the 
occurrence  of  gas  bubbles  can  be  reduced  by  operating  the  sensor  at  a 
very  low  temperature. 

However,  as  the  operating  temperature  of  the  sensor  is  reduced, 
the  sensor  becomes  relatively  more  sensitive  to  temperature 
fluctuations.  Therefore,  the  temperature  of  the  present  water  flow 
facility  will  have  to  be  controlled. 


Thus,  the  use  of  a  hot  film  constant  temperature  anemometer  will 
require  the  construction  of  a  calibration  fac: . ity  which  includes 
temperature  control.  Additionally,  modification  will  have  to  be  made 
to  the  existing  water  tunnel  facility.  These  modifications  include  the 
addition  of  components  to  provide  temperature  control  and  removal  of 
dissolved  gas. 

2.4.2  Proposed  Measurements 

We  propose  to  measure  the  time  resolved  axial  velocity  component 
at  various  radial  positions,  at  the  exit  plane  of  both  a  five  foot  and 
a  ten  foot  tube.  These  measurements  will  provide  a  fundamental  check 
of  the  model  and  also  indicate  with  respect  to  only  the  axial  velocity 
profile  whether  the  flow  is  fully  developed.  We  will  also  develop  the 
instrumentation  techniques  that  will  permit  the  simultaneous  time 
resolved  measurement  of  both  the  axial  and  radial  velocity  components. 


2.5  Heat  Transfer 


jji  Air  Teat  Rig 


One  of  the  objectives  of  this  program  that  is  necessary  to 
obtaining  an  understanding  of  the  fluid  mechanics  in  the  spiral  fluted 
tube  is  the  mesurement  of  the  turbulence  parameters  of  the  flow  field. 

At  the  onset  of  the  program  it  was  thought  that  this  would  best  be 
accomplished  with  hot  wire  anemometry  in  air.  However,  the  subsequent 
results  of  heat  transfer  measurements  at  HTRI  showed  considerable 
en  ancement  with  water  flow  but  not  improvement  with  air.  In  both 
cases  the  direction  of  heat  flow  was  radially  inward.  In  the  HTRI 
experiments,  condensing  steam  at  relatively  high  temperature  was  used 
as  the  heat  source  on  the  outside  of  the  tube,  and  it  was  thought  that 
the  value  of  the  Richardson  Number  was  too  negative.  In  the  Mizushina 
experiments  the  ratio  of  the  turbulent  exchange  coefficients  of  heat 
to  momentum  reached  a  peak  of  3  at  a  slightly  negative  value  of  the 
Richardson  Number  and  decreased  as  this  parameter  decreased.  An 
alternate  explanation  of  the  different  heat  transfer  behavior  in  gas 
as  compared  with  liquid  is  in  the  viscosity  variation  with  temperature. 

In  gases  the  viscosity  increases  with  increasing  temperature,  while  in 
liquids  it  decreases.  It  was  decided  to  examine  this  question  in  more 
detail  by  performing  two  additional  experiments.  The  first  experiment 
allowed  the  heat  flux  to  be  varied  to  include  much  lower  heat  fluxes 
than  were  achieved  in  the  HTRI  experiments.  In  the  second,  the  direction 
of  heat  flow  would  be  reversed  to  radially  outward. 

The  first  experiment  was  completed  during  the  reporting  period  in 
an  air  test  rig  that  was  constructed  and  is  shown  in  Figure  11. 

Filtered  laboratory  air  was  used  as  the  test  fluid  in  an  open  loop 
system.  Steady  air  flow  was  provided  to  the  test  section  by  means  of  a 
centrifugal  air  compressor  in  line  with  a  series  of  flow  straighteners, 
a  contraction  section,  and  a  3-inch  I.D.  aluminum  tube.  Upon  leaving 
the  test  section,  the  fluid  expanded  through  a  sparger  tube  into  a 
second  3  inch  aluminum  tube,  and  was  finally  exhausted  through  a  second 
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contraction  section.  An  inclined  manometer  was  used  to  aeasure  the 
pressure  drop  across  the  first  contraction  to  obtain  the  flow  rate*  and 
a  second  sunometer  Manured  gauge  pressure  in  the  3- inch  tube*  for 
determination  of  fluid  density.  The  momentum  and  continuity  equations 
were  used  to  calculate  fluid  flow  rates.  Test  fluid  temperatures  were 
measured  with  calibrated  platinum  resistance  temperature  probes, 
located  at  the  entrance  to  the  test  section  and  at  the  system  outlet. 

The  outer  surface  of  the  spiral  tube  was  packed  with  heat  conducting 
grease  and  then  installed  inside  a  copper  sleeve*  which  made  complete 
contact  with  the  grease.  Test  section  heating  was  accomplished  by 
wrapping  the  copper  sleeve  connected  to  a  variable  voltage  source 
(0-240  volt)  with  a  resistance  heating  element  (0-1200  watt).  The  heat 
sink,  grease,  and  copper  sleeve  were  used  to  smooth  out  the  heat 
distribution  delivered  by  the  heating  element.  Voltage  across  the 
heater  and  current  through  it  were  measured  for  computation  of  power 
input.  Special  thermal  insulation  covered  the  experimental  system 
from  the  first  contraction  section  to  the  system  exhaust.  Temperature 
distribution  along  the  outer  surface  of  the  spiral  tube  was  measured  by 
means  of  copper-constantan  thermocouples*  with  a  32°  reference, 
installed  in  opposing  sets  of  two  at  a  6-inch  intervals  along  the  tube. 

Data  for  three  different  heat  input  rates  with  3  to  5  flow  rates 
for  each  heat  rate  were  manually  recorded.  These  data  included  inlet 
and  outlet  test  fluid  temperature,  thermocouple  voltage  output,  manometer 
readings  for  pressure  drop  and  pressure  in  the  3-inch  section*  and 
voltage  readings  for  power  input  to  the  heater. 


2.5.1  Data  Reduction 


The  rate  of  heat  input  (Q)  to  the  test  fluid  was  calculated  from 
the  relation 

0-  WC„  (T2-  T{)  (1) 

where  V  (lb/hr.)  is  mass  flowrate  of  the  test  fluid,  Tj(0F)  and  T2(°F) 
are  temperatures  of  the  test  fluid  at  the  test  section  inlet  and  system 
outlet,  respectively,  and  (Btu/lbm°F)  is  the  specific  heat  at 
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constant  pressure  of  the  test  fluid.  The  difference  between  the  rate 
of  beat  input  calculated  in  equation  1  and  that  calculated  from 
voltage- current  measurements  of  the  heating  element  gives  an  estimate 
of  heat  loss  from  the  system  through  pathways  other  than  the  test 
fluid. 


After  an  initial  distance  from  the  test  section  inlet*  skin 
temperature  of  the  tube  increases  linearly  with  longitudinal  distance. 
An  equation  describing  this  linear  relationship  takes  the  form 

t  —  mx  +  b  (2) 

where  t(°F)  =  tube  skin  temperature 
x(ft)  =  distance  along  tube 

A  least  squares  fit  of  temperature  data  obtained  from  the  thermocouples 
was  used  to  determine  the  slope  m  and  the  Intercept  b. 

An  overall  heat  transfer  coefficient  U  (Btu/hr-f t2oF)  can  be 
obtained  from  the  defining  equation 

Q-AUtom  (3) 

where  Q  (Btu/hr.)  =  heat  transfer  rate  (Eq.  1) 

A  (ft^)  =  tube  heat  transfer  surface  area 


Afm  (°F) 


A  tin  Af0U| 

In  (&{;„/ A  foui) 


(4) 


where  At  is  the  log  mean  temperature  difference.  At.  is  the  difference 
si  in 

between  the  test  section  Bkin  temperature  (Eq.  2)  and  the  test  fluid 
bulk  temperature  at  the  test  section  inlet*  and  Atout  is  similarly 
defined  at  the  outlet.  Values  of  0  were  plotted  in  the  form  1/U  versus 
1/V0*8*  the  standard  Vilson  plot  technique*[4]  where  v  is  test  fluid 
velocity  in  the  tube.  A  least  squares  fit  of  the  points  was  obtained 

giving  the  slope  and  intercept  of  the  equation 
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I/U-  mV08  +b 


(5) 


The  value  1/U  Is  actually  a  sum  of  thermal  resistance  due  to  the  tube 
wall  thickness  and  the  fluid  film  thickness  on  the  Inside  wall.  The 
thickness  of  this  film  decreases  with  Increasing  fluid  velocity.  The 
intercept  on  the  vertical  axis  of  the  Wilson  plot,  which  corresponds 
to  infinite  velocity,  represents  the  resistance  of  the  tube  wall  and 
resistances  outside  the  wall  which  in  this  experiment  is  the  leakage 
heat.  A  new  line  can  then  be  constructed  parallel  to  the  original  plot 
but  with  intercept  at  the  origin.  This  line  is  the  resistance  due  to 
the  fluid  film  alone,  and  is  the  inverse  of  the  overall  coefficient 
independent  of  the  other  resistance  terms.  Thus  by  setting  the  inter¬ 
cept  b  equal  to  zero  in  equation  (5) ,  the  resulting  equation  takes  the 
form 


l/hc  =  mv~ot  (6) 

which  is  used  to  obtain  hc«  This  value  is  then  independent  of  thermal 
resistances  other  than  those  caused  by  the  structure  of  fluid  flow 
inside  the  tube.  Further,  by  taking  a  least  squares  fit  in  constructing 
the  original  plot,  correction  is  made  for  the  slight  scatter  of  data 
resulting  from  variations  in  the  rate  of  heat  input  as  the  flow  rate 
is  varied. 

The  Nusselt,  Reynolds  and  Frandtl  numbers  are  dimensionless 
quantities  defined  by 

Nu  -  h  D. /k 
n 

Re  -  v  D. /v 

n 

Pr  ■  c  y/k 
P 

where  k,  v  and  y  are  evaluated  at  mean  bulk  fluid  temperatures. 
Experimental  results  were  plotted  in  the  form  Nu*  versus  Re  where 

Nu *- Nu/[Pr°<  (  n  /  Fj° 5  ]  . 
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T,  is  the  arithmetic  mean  bulk  fluid  temperature,  and  T  is  the 

D  V 

mean  wall  temperature.  Finally,  a  correlation  to  the  heat  transfer 
performance  of  a  smooth  tube  was  given  by  the  well  known  Dittus  Boelter 
equation [5] 

Nu*  -  0.023  Re0  8 

The  results  are  shown  in  Figure  12  where  the  Nusselt  Modulus  vs. 
Reynolds  Number  variation  is  shown  for  the  three  heating  rates  and  the 
Dittus  Boelter  correlation  is  shown  as  a  straight  line.  The  data  agree 
with  the  straight  tube  correlation  with  the  HTRI  data  at  the  higher 
heating  rates.  Table  1  gives  the  tabulated  values  of  the  measurements 
and  calculations. 
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HEAT  TRANSFER  AIR  TEST 


3 .  THEORETICAL  PROGRAM 


3.1  lau.flduat.ittn 

The  spirally-fluted  tube  appears  to  have  attractive  features  relative 
to  enhanced  heat  transfer  and  thus  to  heat  exchangers.  In  particular*  the 
experimental  data  indicate  increased  heat  transfer  without  a  corresponding 
increase  in  pressure  drop.  From  a  fundamental  point  of  view,  these  features 
call  for  an  explanation  of  the  mechanism  by  which  the  heat  transfer  at  the 
wall  is  facilitated  while  the  momentum  transfer  leading  to  pressure  drop  is 
not  similarly  increased.  It  is  worth  noting  that  analogous  features  are 
observed  in  geophysical  turbulence  where  gravity  provides  a  body  force  field 
corresponding  to  the  centrifugal  force  field  associated  with  swirl 
(cf..e.g.*  Meroney  [7]).*  It  is  the  purpose  of  this  report  to  outline  our 
current  approach  to  the  analysis  of  the  flow  within  spirally-fluted  tubing. 

It  will  become  readily  evident  as  our  discussion  proceeds  that  the 
analysis  of  turbulent  flow  in  such  tubing  cannot  be  fruitfully  developed 
without  an  accompanying  experimental  effort  providing  detailed*  accurate  and 
extensive  data  on  the  statistics  of  the  velocity  components  and  subsequently 
of  the  temperature  and  the  temperature  velocity.  Global  measurements  of 
heat  transfer  and  pressure  drop*  the  only  measurements  carried  out  to  date 
in  spirally-fluted  tubing*  are  not  sufficient  to  guide  an  analysis  whose 
purpose  is  to  provide  more  than  correlations  of  experimental  results*  namely 
a  fluid  mechanical  explanation  of  the  transfer  mechanisms  involved.  Such  an 
explanation  is  necessary  for  the  optimal  design  of  such  tubing. 


*Meroney  makes  the  following  informative  statement:  "in  the  unstable  case* 
a  large  flat  heat  flux  is  established  quite  rapidly*  momentum  transport  may 
be  significantly  smaller  in  proportion....  It  appears  that  buoyancy¬ 
generated  eddies  cause  relatively  little  momentum  transport*  but  they  are 
quite  effective  at  carrying  thermal  energy." 


The  need  for  experimental  data  is  dramatically  clear  when  the  boundary 
conditions  to  be  applied  at  the  surface  of  the  flutes  are  confronted.  The 
flow  within  and  immediately  adjacent  to  the  flutes  must  play  an  important 
role  in  determining  tube  performance.  It  is  very  complex  and  is  being 
studied  separately  in  the  continuing  work.  Under  these  circumstances*  such 
a  role  must  be  reflected  in  a  modification  of  the  “law  of  the  wall*~  a  modi¬ 
fication  which  can  only  be  determined  by  appropriate  experiment.  The 
possible  importance  of  the  wall  region  is  indicated  by  the  provisionally 
observed  difference  between  tube  performance  with  water  versus  air  as  the 
working  fluid;  the  performance  with  the  former  appears  to  be  significantly 
better  than  with  the  latter  even  within  the  same  range  of  the  parameters 
usually  considered  to  determine  flow  behavior. 

In  an  earlier  study  [8].  we  apply  the  results  of  So  [9]  to  describe  the 
fully-developed  turbulent  flow  in  tubing  with  spiral  flutes.  So  assumes 
that  in  axisymmetric  swirling  flow*  the  various  production  and  dissipative 
terms  in  the  describing  equations  are  so  dominant  that  the  equations  are 
essentially  satisfied  by  equating  these  terms.  The  results  involve 
algebraic  expressions  for  the  turbulent  exchange  coefficients  in  the  stream- 
wise  and  radial  directions.  To  utilize  these  relations*  it  is  necessary  to 
assume  the  radial  distribution  of  a  turbulent  length  scale.  So  argues  that 
this  scale  is  unaltered  by  swirl;  thus  in  [3]  we  assume  the  standard  length 
scale  distribution  l(r)  taken  from  [10].  In  [11],  So  extends  this  approach 
to  flows  with  heat  transfer. 

It  is  widely  accepted  that  Reynolds  stress  closure  represents  the  more 
complete  formulation  for  turbulent  shear  flows  and  that  as  this  closu/. 
methodology  evolves*  it  will  result  in  a  standard  procedure  for  engineering 
calculations  applicable  to  such  flows.  Accordingly*  in  contemplating  an 
approach  for  calculating  not  only  the  velocity  field  in  spirally-fluted 
tubing  but  the  thermal  field  as  well*  it  appears  that  the  full*  second 
moment  treatment  of  turbulent  flows  will  provide  a  more  general  framework 
for  the  incorporation  of  experimental  data  than  the  analysis  of  So.  The 
basis  for  such  a  treatment  is  provided  by  Launder  and  Morse  [12]  who  apply 


Reynolds  stress  closure  to  axisymaetric  shear  flows  without  and  with  swirl. 
Although  it  is  possible  to  undertake  a  detailed  assessment  of  the  extensive 
modeling  called  for  in  their  treatment*  it  is  consistent  with  our  view  on 
the  need  for  an  evolutionary  approach  to  accept  their  formulation  as  a 
starting  point.  In  this  regard*  it  should  be  noted  at  the  outset  that  Ref. 
[7J  is  in  the  nature  of  a  progress  report;  it  indicates  that  in  its  present 
stage  of  development*  the  theory  of  axisynmetric  shear  flows  is  not 
adequate.  Launder  and  Morse  show  that  the  spreading  rate  for  nonswirling 
jets  is  predicted  to  be  50Z  too  high.  In  addition*  with  swirl  the  shear 
stress  .  .  uw  has  the  wrong  sign  and  thus  the  numerical  solutions  display 
a  reduced  rate  of  spread  in  contrast  to  the  strong  augmentation  found  in 
practice."  These  shortcomings  suggest  the  difficulties  of  predicting 
swirling  turbulent  flows  and  the  importance  of  close  coordination  between 
experimental  and  predictive  efforts.  We  thus  adopt  the  view  that  despite 
current  deficiencies,  the  theory  can  in  due  course  evolve  to  provide  a 
useful  predictive  capability. 

As  a  further  justification  for  our  adopted  approach*  we  note  that 
second -moment  closure  exposes  in  a  clear  fashion  the  mechanism  whereby  swirl 
affects  the  heat  transfer  characteristics  of  the  tube.  We  find  that  the 
term  0'2dp/dr  *  the  product  of  the  intensity  of  the  temperature  fluctuations 
and  the  radial  pressure  gradient  describes  that  mechanism.  It  is  interest¬ 
ing  to  observe  that  this  term  is  the  counterpart  of  the  gravity  term  in  the 
atmospheric  boundary  layer. 

At  present*  attention  focuses  on  analyzing  the  velocity  field  within 
the  tube.  The  coupling  with  the  temperature  field  is  considered  as  a  second 
step.  This  sequential  approach  implies  that  we  deal  with  small  rates  of 
heating  and  cooling. 


r 
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3.2  Aaalyaia 

Figure  13  taken  from  C 83  shows  schematically  the  flow  and  coordinate 
system  considered.  We  assume  the  flow  is  of  constant  density  and  fully- 
developed.  i.e..  only  the  mean  pressure  "p  varies  in  the  streamwise 
direction.  We  follow  the  practice  in  [8]  and  [12]  and  denote  the  mean 
velocity  components  by  capital  letters  and  the  related  fluctuations  by  lower 
case  letters.  Thus,  the  mean  velocity  components  of  interest  are  U(r)  and 
W(r)  with  fluctuations  u  and  w.  Although  the  mean  radial  velocity  is  zero 
within  the  core  of  the  tube,  its  fluctuation  v  plays  an  important  role. 

3.2.1  Describing  Equations 

With  these  assumptions,  preliminary  work  indicates  that  we  must 
consider  the  following  equations  resulting  directly  from  [12]: 
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Conservation  of 
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In  these  equations,  the  notation  is  standard  but  several  quantities  call  for 
comment.  The  quantity  k  is  the  turbulent  kinetic  energy  l/2u^u^  while  e  is 
the  velocity  dissipation.  The  units  of  k/e  are  (L/U)  so  that  the  coeffic¬ 
ients  in  the  diffusion  terms  have  the  units  of  (UL)  as  required. 

We  need  additional  equations  to  complete  the  formulation,  but  it  is 
convenient  to  discuss  Eqs.  (l)-(5)  before  proceeding.  In  these  equations, 
the  following  empirical  constants  appear:  c^;  C2  with  the  related  constants 

a  =  (8  +  c;>)/l  1,  /3  =*  (8c2  —  2)/ll  and  y  —  (30c2  -  2)/55;  cs\  re);  and  re¬ 
values  of  these  constants  are  given  in  C12]. 

In  the  absence  of  swirl.  W  =  0  and  Eqs.  (l)-(5)  represent  five 
equations' for  the  three  velocity  intensities,  u^.  v^»  w^,  the  mean  velocity 
U  and  the  dissipation  c.  Solutions  to  these  equations  subject  to  regularity 
conditions  at  the  axis  (r  =  0)  and  boundary  conditions  at  the  wall  describe 
the  turbulent  flow  in  either  unfluted  tubes  or  tubes  with  straight  fluteB. 
Clearly  such  solutions  must  be  obtained  as  a  first  step  in  the  development 
of  our  analysis. 

The  continuity  and  three  mean  momentum  equations  do  not  enter 
explicitly  in  our  formulation;  the  first  of  these  equations  leads  to  V  =  0 
while  the  x-wise  momentum  equation  becomes 
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where  A  =  2uv(r  =  R)/U^  is  the  friction  factor  and  0  is  the  mean  velocity 
defined  by 
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m 


(6) 


He  thus  obtain  the  radial  distribution  of  Reynolds  shear  stress 


u\  - 


(7) 


so  that  uv  is  not  a  dependent  variable  in  Eqs.  (l)-(5).  The  radial  momentum 
equation  determines  the  radial  pressure  gradient  after  the  characteristics 
of  the  velocity  are  known;  we  have 
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As  suggested  earlier,  Eq.  (8)  will  enter  significantly  in  the  analysis  of 
the  thermal  field  in  the  tube  but  for  the  present  may  be  set  aside. 

Finally,  the  aximuthal  momentum  equation  leads  to  the  important  result 
that  the  Reynolds  shear  stress  vw  =  0.  He  thus  find  that  only  two  of  the 
possible  three  shear  stresses  are  operative  in  the  flow  under  consideration. 

Uith  swirl  one,  additional  equation  is  required,  that  for: 
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We  must  now  take  up  the  issue  of  the  radial  distribution  of  W;  physical 
intuition  suggests  that  within  the  central  tube*  there  is  a  central  core 
having  a  solid  body  rotation  so  that  W  =  ftr  where  ft  is  a  measure  of  the 
extent  of  rotation.  This  core  region  extends  to  the  neighborhood  of  the 
flutes  in  the  case  of  fluted  tubing.  Thus*  the  radius  r  =  R  takes  on  two 
meanings.  For  a  smooth  tube*  R  represents  the  interior  radius.  Such  a  tube 
may  be  fixed  so  that  W  =  0  or  rotating  with  an  angular  velocity  ft  .  Experi¬ 
mental  results  relative  to  the  mean  velocity  distributions  in  the  latter 
situation  [13]  and  [14]  indicate  results  similar  to  those  predicted  in  [8]. 
namely*  profiles  consistent  with  the  stabilizing  effect  of  rotation.  The 
absence  of  turbulence  measurements  vitiates  the  value  of  these  results. 

For  fluted  tubes*  the  radius  R  should  be  considered  R  -  nh  as  indicated 
in  Figure  13  and  is  considered  to  determine  the  edge  of  the  core  involving 
rigid  body  rotation.  The  influence  of  the  flutes  on  the  flow  is  reflected 
in  the  boundary  conditions  to  be  imposed  at  r  =  R.  The  proper  formulation 
of  these  conditions  must  be  guided  by  experiment.  With  W(r)  assumed  given. 
Eq.  (9)  completes  the  system  of  six  equations  in  six  unknowns. 


3.2.2  Boundary  Conditions 

The  boundary  conditions  at  the  axis  r  =  0  are  given  by  the  physical 
requirements  of  regular  behavior.  We  thus  have: 

i.o 

dr  dr  dr  dr  dr 

U  -  U0,  uw  -  0  (10) 


In  Eqs.  (l)-(5)  and  (9)*  there  are  three  intensities*  the  dissipation  and 
the  uw-stress*  each  involving  second  derivatives*  and  the  U-velocity 
component  involving  a  first  derivative  for  a  total  order  of  eleven.  The 
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condition  on  (dU/dr)  is  automatically  antiaified  if  the  other  aix  conditions 
are  applied  ao  that  there  are  four  free  quantities  at  r  =  0  to  be  selected 
so  that  four  boundary  conditions  at  r  =  R  can  be  imposed.  We  find  it 
convenient  to  let  u^(0)»  v^(0)  and  e(0)  be  three  of  the  quantities  in 
question  with  w^(0)  =  v^(0)  the  fourth. 

In  contrast*  the  boundary  conditions  to  be  imposed  at  the  edge  of  the 
core  in  order  to  reflect  the  influence  of  the  flutes  on  the  flow  are 
uncertain.  Before  making  even  a  provisional  suggestion  concerning  the 
nature  of  these  conditions*  it  is  worth  discussing  several  possible  notions 
of  the  flow  behavior  in  the  neighborhood  of  the  flutes.  Suppose  the  resul¬ 
tant  mean  velocity  due  to  the  combination  of  the  U-  and  W-velocity  compon¬ 
ents  is  aligned  with  the  flutes  so  that  for  r  ~  R  ~  W  ~  U  tan  i.  Within  the 
immediate  vicinity  of  the  flutes*  both  W  and  U  decrease  and  approach  aero. 
The  distribution  of  azimuthal  velocity  within  the  flutes  is  thus  destabil¬ 
izing  so  that  we  think  of  the  flute  region  as  a  region  of  high  turbulence 
intensity  tending  to  counteract  the  stabilizing  influence  of  swirl  in  the 
core  region. 

If  the  resultant  mean  velocity  in  the  neighborhood  of  the  flutes  is  not 
aligned  with  the  flutes*  the  flow  will  lead  to  further  increases  in  turbu¬ 
lent  intensities.  This  discussion  suggests  one  obvious  measurement*  namely 
that  of  the  radial  distribution  of  the  two  mean  velocities*  U  and  W*  with 
special  attention  devoted  to  the  neighborhood  of  the  flutes. 

If  we  adopt  the  view  that  prediction  of  tube  flow  without  swirl  is  the 
first  order  of  business*  the  boundary  conditions  at  r  *  B  are  clearly  given 
by  the  law  of  the  wall.  In  fact*  this  law  provides  a  starting  point  for  a 
provisional  formulation  which  can  be  applied  to  fluted  tubes  pending  experi¬ 
mental  data  leading  to  an  improved  description  of  conditions  in  the  neigh¬ 
borhood  of  the  flutes.  Thus*  we  require  as  r  ~  R 


U-{m)^[A  «((*  -r)  uvm/v\  +  B  . 


42 


By  eliminating  the  distance  from  the  core  radius  R  -  r,  this  equation  can  be 
rewritten  as 


U-(i Tv) 1/2 
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uvA 

v(-dU/dr) 


0  . 


(11) 


Equation  (11)  assures  that  the  velocity  and  the  velocity  gradient  in  the 
radial  direction  at  the  edge  of  the  core  are  consistent  with  the  log  portion 
of  the  law  of  the  wall.  He  expect  that  this  equation  may  prevail  even  with 
swirl  but  are  prepared  to  make  modifications  if  experimental  results  so 
indicate. 


The  other  two  conditions  to  be  imposed  at  r  =  R  are  more  uncertain;  the 
work  of  Launder  et  al.  [15]  suggests  four  possible  choices,  three  relating 
u^,  v^»  and  w^  to  A  and  a  fourth  relating  e  to  A  dU/dr.  In  the  absence  of 
guiding  experimental  data,  we  impose  a  condition  on  the  intensity  of  the 
radial  velocity  fluctuations  v^  since  intuition  suggests  that  the  most 
direct  influence  of  the  flutes  on  the  core  flow  is  likely  to  be  via  that 
intensity.  Thus,  at  r  =  R.  we  take 

v2-K2«v  (12) 


where  *n  empirical  constant.  As  the  final  condition,  we  adopt  the 

condition  on  the  dissipation  from  [15].  namely 


(  -  - 


(13) 


Although  we  have  emphasised  the  provisional  nature  of  the  three 
boundary  conditions  at  r  *  t,  it  should  be  noted  that  the  general  structure 
of  the  analysis,  i.e..  the  availability  of  three  quantities  at  r  =  0  so  that 
three  conditions  at  r  s  R  can  be  specified,  is  firm.  Accordingly,  other 
conditions  can  be  imposed  without  altering,  in  a  significant  fashion,  our 
analysis. 
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With  swirl*  we  need  another  condition  «t  r  =  R.  If  the  ween  flow  at 


the  edge  of  the  core  is  taken  to  be  tangent  to  the  flutes*  we  have 

W  -  U  tan  *•  HR  .  (14) 

In  connection  with  Eq.  (14)*  we  note  that  nuaerical  analysis  is  facilitated 
if  an  inverse  method  is  followed*  thus  we  specify  ft  and  calculate 
a  posteriori  the  appropriate  flute  angle 

3.2.3  Hondimensionalization  and  Reduction  to  First  Order  Equations 

We  use  the  centerline  mean  velocity  U0  and  the  core  radius  R  as  the  two 
nondimensionalizing  quantities.  Accordingly*  we  define  the  following  depen¬ 
dent  variables: 

MU-VIU0,  Iu  -»  u*/U$  ,  h-S/Ui,  Iw  -  w2/U^ 

,  ,  (15) 

E  -  tR/ U$  ,  Fuw  -  uw/Ut ,  K  -  y  (/„  +  /v  +  /„) 


In  order  to  reduce  the  differential  equations  to  first  order*  we  introduce 
the  following  additional  dependent  variables: 
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where  prime  denotes  differentiation  with  respect  to  n  =  r/R*  the  single 
independent  variable. 

With  the  new  variables*  Sqs.  (l)-(5)  and  (9)  become: 
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where  T  =  El  / nK  is  a  dimensionless  transfer  coefficient  and  £1  £  tJR/U  is 
v  o 

a  dimensionless  swirl  parameter.  These  equations  are  to  be  solved  subject 
to  boundary  conditions  at  i)  *  0: 
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where  NR  s  U  R/v  is  the  usual  Reynolds  number  based  on  the  mean  velocity. 
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The  solutions  to  these  equations  depend  on  parameters  X ,  Ng  and  ft  and 
on  various  empirical  constants  which  we  can  assume  to  be  known  but  subject 
to  adjustment  to  bring  prediction  and  experiment  into  agreement.  The  speci¬ 
fication  of  a  value  for  X  calls  for  further  empiricism  relating  the  pressure 
drop  characteristics  to  the  Reynolds  number*  e.g.»  to  UR/v,  and  to  the  flute 
geometry*  i.e.,  to  the  height*  pitch*  and  spiral  angle  of  the  flutes.  In 
[3]*  we  use  the  pressure  drop  equation  due  to  filasius  written  in  terms  of 
the  hydraulic  radius  and  relate  that  radius  to  the  flute  geometry.  The  same 
treatment  can  be  applied  to  the  present  analysis  so  that  we  compute  as  part 
of  each  solution  the  ratio  U/U0  using  Eq.  (6). 

3.2.4  Strategy  of  Solution 

The  numerical  solution  of  this  imposing  set  of  eleven  differential 
equations  involving  eight  boundary  conditions  at  n  =  0  and  three  at  n  =  1 
clearly  involves  two-point  boundary  values.  He  propose  to  use  the  method  of 
quasilinearization  to  obtain  the  solutions.  In  this  approach  integration  of 
the  full*  nonlinear  equations  with  assumed  values  of  Iu(0)>  Iy(0)  and  E(0) 
plus  three  homogeneous  equations  obtained  by  quasilinearization  is  initiated 
at  n  =  0  and  carried  out  to  o  =  1.  The  arbitrary  constants  which  can 
multiply  the  solutions  to  the  homogeneous  equations  are  selected  so  that  the 
boundary  conditions  at  n  =  1  are  satisfied  and  so  that  the  values  of  the 
guessed  quantities  are  systematically  improved.  Convergence  is  achieved 
when  no  further  improvements  are  called  for. 


3 .3  Caac India* . teuiln 

We  have  outlined  an  analysis  applicable  to  the  fully  developed 
turbulent  flow  in  tubes  with  spiral  fluting.  The  general  approach  involves 
a  second  moment  closure  which  is  believed  to  provide  a  formulation  permit¬ 
ting  incorporation  of  essential  experimental  results.  The  importance  of 
such  results  in  the  development  of  boundary  conditions  to  be  applied  in  the 
neighborhood  of  flutes  is  emphasized.  The  present  analysis  utilizes*  in  a 


highly  tentative  way*  conditions  which  complete  the  mathenatical  formulation 
and  which  have  some  physical  content.  However*  we  expect  that  these  condi¬ 
tions  will  require  modification  when  appropriate  experimental  results  are 
available.  Such  modification  can  be  incorporated  in  our  general  approach 
without  difficulty. 


4.  CONTINUING  WORK 


The  work  on  the  experimental  and  analytical  programs  on  the  Fluid 
Mechanics  and  Heat  Transfer  of  Spirally  Fluted  Tubing  is  continuing  and 
will  be  reported  on  in  subsequent  reports.  The  principle  areas  are 
listed  below. 


4.1  Kxpprimpntal  Program 

4.1.1  Flow  Visualization  Studies 

The  flow  visualization  studies  are  being  directed  toward 
obtaining  quantitative  data.  As  was  mentioned  earlier*  this  requires 
flash  synchronization  with  a  high  speed  motion  picture  camera.  The 
elements  of  the  system  have  been  breadboarded  and  are  awaiting  the 
installation  of  the  new  cast  section  in  the  Water  Tunnel. 

4.1.2  Turbulent  Measurements 

Hot  film  measurements  are  planned  in  the  water  tunnel  after 
completion  of  the  flow  visualization  studies.  Work  is  in  progress  on 
the  calibration  rig. 

4.1.3  Heat  Transfer  Measurements  with  Air 

A  new  air  test  rig  has  been  designed  and  the  parts  manufactured 
for  the  measurement  of  the  overall  heat  transfer  coefficient  with  the 
direction  of  heat  radially  outward.  The  equipment  is  being  assembled 
and  the  tests  will  be  performed  shortly. 
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4.2  Theoretical  Program 

4.2.1  Continuing  Hfllk 

The  analysis  described  in  Section  3.0  is  continuing  and  the 
programming  of  the  equations  Is  in  progress. 

4.2.2  Additional  Theoretical  Program 

A  subcontract  has  been  issued  to  Prof.  B.  E.  Launder  at  the 
University  of  Manchester  to  develop  the  methodology  for  the  calculation 
of  the  flow  field  in  and  adjacent  to  the  flutes.  Work  has  started  with 
contributions  from  Dr.  A.  D.  Gosman  of  the  University  of  London  and  Dr. 
Bergeles  of  the  University  of  Athens. 
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LIST  OF  SYMBOLS 

Q  Rate  of  heat  flow  (Btu/hr.) 

W  Maas  flow  rate  (lbm/hr.) 

Cp  Specific  heat  at  constant  pressure  (Btu/lbmOF) 
t  Spiral  tube  outer  wall  temperature  (°F) 

U  Overall  heat  transfer  coefficient  (Btu/hr-f t^°F) 

Uc  Corrected  overall  coefficient  ( Btu/hr- ft^F) 

Dn  Hydraulic  diameter  (ft.) 

K  Test  fluid  thermal  conductivity 

v  Kinematic  viscosity  (ft^/hr.) 

u  Absolute  viscosity  (lb/hr-ft.) 

V  Average  test  fluid  velocity  (ft. /hr.) 
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